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Abstract: We use the mapping between Burgers' equation and the prob- 
lem of a directed polymer in a random medium in order to study the fully 
developped turbulence in the N dimensional forced Burgers' equation. The 
stirring force corresponds to a quenched (spatio temporal) random potential 
for the polymer. The properties of the inertial regime are deduced from a 
study of the directed polymer on length scales smaller than the correlation 
length of the potential, which is not the regime usually considered in the case 
of polymers. From this study we propose an Ansatz for the velocity field in 
the large Reynolds number limit of the forced Burgers' equation in N dimen- 
sions, which should become exact in the limit N oo. This Ansatz allows 
us to compute exactly the full probability distribution of the velocity differ- 
ence u{r) between points separated by a distance r much smaller than the 
correlation length of the forcing. We find that the moments < u'^{r) > scale 
as r'^^"^ with ((q) = 1 for all g > 1 (in particular the g = 3 moment agrees 
with Kolmogorov's scaling C(3) = 1). This strong 'intermittency' is related 
to the large scale singularities of the velocity field, which is concentrated on 

-'^Unite propre du CNRS, associee a I'Ecole Normale Superieure et a I'Universite de 
Paris Sud 



a iV — 1 dimensional froth-like structure, which is turn related to the one- 
step replica symmetry broken nature of the associated disordered problem. 
We also discuss the similarities and differences between Burgers turbulence 
and hydrodynamical turbulence, and we comment on the anomalous tracer 
fluctuations in a Burgers turbulent field. This replica approach being rather 
unusual in turbulence problems, we provide all the necessary details of the 
method. 
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1 Introduction 



The theory of statistical turbulence is more than fifty years old, but its status 
is still not fully satisfactory. The simplest and surprisingly robust approach 
is Kolmogorov's dimensional analysis, which leads to the celebrated k~^^^ law 
for the velocity cascade 0]. However, analytical calculations immediately 
lead to difficulties: the simplest closure scheme to deal with the non linear 
term in the Navier-Stokes equation ('DIA') does not reproduce Kolmogorov's 
scaling [^, ^, ^. A lot of effort has been devoted, in particular by Kraichnan, 
to understand why this was so. More refined schemes were proposed to 
recover the /c~^/^ law and even to calculate adimensional constants p, ||, 0], 
but they are based on uncontrolled assumptions (like the RG approach), or 
lead to enormous calculations, which are difficult to manipulate and are not 
fully transparent from a physical point of view 

More recently, a tremendous activity has developed on the intermittency 
(or multifractal) corrections to Kolmogorov's scaling: higher moments of 
the velocity field do not seem to scale with the one originally predicted by 
Kolmogorov. Many interesting suggestions have been put forward to describe 
and explain this feature, starting by Kolomgorov himself 0, ||, |^, [1^, |ll], |12 
|T3|,|T 



In recent years, non linear partial differential equations with noise 
have been the focus of quite a number of studies in the context of grow- 
ing interfaces, with the 'KPZ' equation standing out as a new paradigm 
[p!5| , 0, 0. This equation is in fact a variant of Burgers' equation, and 
has (together with its many siblings) a wide range of applications in different 
physical contexts [0. Interestingly, this field is also related (through an 
appropriate 'Cole-Hopf mapping) to the physics of disordered systems, in 
particular elastic strings in random media ('directed polymers'). The Burg- 
ers problem in N dimensions is then equivalent to the problem of a directed 
polymer in + 1 dimensions, which is the space time of the original problem. 
The forcing term in the Burgers equation translates into a random potential 
for the directed polymer, which is quenched in space time: to each realiza- 
tion of the stirring force in the Burgers language corresponds one sample of 
a directed polymer. This enables one to adapt techniques originally devised 
for with spin-glasses and obtain original results on the underlying non-linear 



equation (or vice- versa) |T^, ^ ^ ^ It is the aim of the present 



paper to exploit in detail this mapping, which allows us to propose an orig- 
inal Ansatz for the velocity field for the randomly stirred Burgers equation, 
which should become exact solution in high dimensions and in the limit of 
large Reynolds number. Scaling in the inertial range can then be precisely dis- 
cussed: we obtain in closed form the full probability distribution for velocity 
differences. In particular, the third moment of the velocity difference grows 
linearly with distance, i.e d la Kolmogorov. In fact, Kolmogorov's dimen- 
sional analysis should be directly applicable to Burgers' turbulence. We find 
however very strong intermittency effects, which we relate in a quantitative 
manner to the existence of large scale structures, in the form of singularities 
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concentrated on time dependent hypersurfaces. Although not exact in finite 
dimension, we expect that our (variational) description remains qualitatively 
correct even in one dimension. We also obtain the dynamical evolution of 
the velocity field: we find that the field is convected away by the largest 
structures, corresponding to a dynamical exponent equal to 2; = 1 rather 
than the one obtained from Kolmogorov's scaling, = 2/3. 

In the bulk of the paper, we primarily focus on the 'dictionary' between 
Burgers' turbulence and disordered systems and discuss a number of physical 
points, relegating more technical points to various appendices. We show that 
Kolmogorov's scaling has a counterpart in the 'directed polymer' language, 
where it is known as the Larkin-Ovchinnikov scaling. It is however easy to 
see that this 'naive' scaling cannot hold, as is confirmed by the full calcu- 
lation. We describe in physical terms the nature of the velocity field, and 
argue that there should be a large distance regime (beyond the 'injection' 
scale) characterized by a non-trivial exponent. We briefiy discuss the prob- 
lem of a passive scalar in such a velocity field. In the conclusion, we compare 
our results to other approximation schemes, and comment on the possible 
differences with 'true' Navier-Stokes turbulence. 



2 From Burgers turbulence to directed poly- 
mers: dictionary and dimensional analysis 

The problem we shall consider is that of a randomly forced potential flow in 
N dimensions, which follows the Burgers equation 0: 

- S/)v = uV^v + f{x,t) (2.1) 



dt 



+ {v 



where v is minus the gradient of a velocity potential h{x,t), and / is a 
randomly fluctuating force, which is also minus the gradient of a potential 
(j){x,t). V and x are A^- dimensional vectors, and there is no constraint on 
V.v. Here, we wish to describe the problem of a fluid which is randomly 
stirred only at very large lengthscales. We shall thus take to be Gaussian, 
with fluctuations given by: 



(f , t)0(f' , f) = eA^N6{t - t') exp • 



\x 



X 



'/\2 



2A^A2 



(2.2) 



Here A is the 'injection' length, over which the forcing is roughly constant. 
The correlation time of the forcing was set to zero; a flnite (small) correlation 
time would however not affect the following conclusions. The stirring force 
correlation thus reads: 



Mx,T)r{x',T')=eS{T 



T 



(X — X 



x'r 



iVA2 





(f — x'Y 


exp — 


2NA^ 



(2.3) 



^The density of the fluid is taken equal to one. 
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As will be clear below, this ensures that the injected energy per unit time, 
defined as the increase of kinetic energy is equal to (Note that the 
dimension of e is The dependance of the force correlations on the 

dimension N has been chosen such as to insure the existence of the large N 
limit. 

The typical velocity at the injection scale is, from dimensional consider- 
ations, va = (eA)3, which allows us to define the Reynolds number as: 

ife . ^ = . (2.4) 

We shall be interested in studying the statistics of the velocity field at large 
Re. 

Let us now use the standard techniques to transform this problem into 
a directed polymer. Integrating once Eq. (1), with v — —Vh one finds the 
so-called Kardar-Parisi- Zhang (KPZ) equation: 

dh{x,t) 1,^, ^vv2 , V72 



dt 2 



(V/i(f , t)y + vV'h{x, t) + 0(f , t) (2.5) 



which describes, in particular, the growth of a surface under a random 'rain' 
of particles, the flux of which is given by (f){x,t). At this stage, the crucial 
difference with previous work on the KPZ equation is in the correlations of the 
noise (j){x,t). Wc shall thus be mostly concerned, in the following, with the 
velocity fleld statistics at length scales sma//er than A; for length scales larger 
than A, the KPZ scaling prevails. Now, as is well known, the KPZ equation 
(2.5) can be transformed (through a 'Hopf-Cole' transformation) into a linear 
problem describing a directed line (polymer) in a random potential. Setting 
h{x,t) — 2i'\og Z{x,t), one finds that Z{x,t) obeys the equation 

= uV'Z{x, t) + ^^{x, t)Z{x, t) (2.6) 

which is the equation for the partition function of an clastic string in a 
random potential V{x,t) = ^(f){x,t), subject to the constraint that its end 
point is fixed at {x,t). Said differently, the solution of Eq. (2.3) can be 
written as a path integral (sum over all configurations): 



p{xo)dxo / d[x{T)] exp(-7i) . (2.7) 

Jg(t=0)=xo 

The hamiltonian 7i is given by: 



rt 

n = I (It 





+v{x{t),t) 



{2.1 



The 'dictionary' between the two problems is the following: The temper- 
ature scale of the polymer has been chosen equal to one. Then the elastic 
modulus of the polymer is 

c.l. (2.9) 



4 



The random potential seen by the polymer V has a gaussian distribution, 
with a second moment given by: 



V{x, r)V{x', r') = WNS{t - t') exp 



2iVA2 



(2.10) 



where the strength of potential fluctuations, W, is related to the energy 
density e through: 

W^-^, (2.11) 

and the length scale of fluctuations of the potential is equal to A, the length- 
scale at which the Burgers fluid is stirred. The probability distribution p{xo) 
of the initial point of the polymer is related to the initial conditions of the 
velocity through: 

v{x, t = 0) = -2vVp{x) . (2.12) 

The non-linearity has disappeared from this formulation, and has been 
replaced by the famous problem of disordered systems, which is to average 
the logarithm of a partition function (in order to calculate various moments 
of the velocity field v). In the next section we shall deal with this problem 
using the replica trick and a variational method which becomes exact in the 
limit of very large dimensions N . 

Before turning to this calculation, it is useful to give some kind of qual- 
itative (more or less dimensional) analysis of the directed polymer, and to 
enrich our 'dictionary' by stating its counterpart in the turbulence language. 
The study of an elastic structure like the directed polymer in presence of a 
random potential has been discussed in several works recently [^,21]. In 



the phase where the disorder is strong, one expects a scaling behaviour of 
the lateral fluctuations of the polymer described by a wandering exponent C,: 



< {x{t) - x{t')Y > ~ A|t - tf (2.13) 

The thermal fluctuations are irrelevant at large distance and this scaling 
also holds in the zero temperature limit, or for the disconnected correla- 
tion (< x{t) — x{t') >)^. As for the free energy differences for two polymers 
finishing at point t, x and t', x', they scale as: 

{h{x,t)-h{x',t')f^\x-x'\'^'^gn • (2-14) 

The t ^ t' limit then implies that the free energy difference at points x and 
x' scales as \x — x'p'^/'', and therefore the difference of velocities in Burgers 
equation should scale as: 

{v{x,t)-v{x',t')f ^ \m - mr-'''9. (^y^^) • (2-15) 

If this last scaling form holds, the galilean invariance of Burgers equation 
implies that the two terms in dv/dt and {v ■ V) v should scale in the same 
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way under a rescaling x ^ bx and t b^^^'t, which imphes a; = 2^ — 1 
[0, (Another way to argue about this identity directly on the polymer 
problem is by observing that the fluctuations of elastic energy in the directed 
polymer scale as x So we are left with only one scaling 

exponent. 



An important point which was discussed in [^6], ^ is the existence of two 
distinct scaling regimes, in the case where the random potential has a large 
correlation length, A. The regime which is most studied for directed polymer 
is the large time regime where the transverse fluctuations of the polymer are 
much larger than A. But there also exists a short distance regime where the 
transverse fluctuations of the polymer are smaller than the correlation length 
of the potential. This is obviously the regime which interests us most for the 
turbulence problem (although we shall return to the long distance regime 
later on). This regime holds for time differences shorter than a typical time 
scale r*, defined from: 



< {x{t + T*) - x{t)f > - A"^ . (2.16) 

One expects to be allowed to linearize the random potential in this slowly 
varying regime, and to study the much simpler (linear) problem of a polymer 
with a random force, defined by the Hamiltonian: 



rt 

n = I dr 







(2.17) 



2 \dT^ 

(where the random force is of order If^l ~ ^/W/A). This random force 



problem was studied by Larkin and Ovchinikov long time ago [^. In this 
regime there is no metastable state and the problem is easily solved for one 
given sample. Assuming for simplicity periodic boundary conditions, one 
obtains the Fourier transform of the average polymer's position:: 

<f(^)>=M. (2.18) 

C 00^ 

A dimensional analysis of the Fourier transform would lead to: 

|<f(t)-f(0>P^^|t-tf, (2.19) 

from which one deduces the scaling exponent in this regime, and the value 
of T*: 

C, = -, ., = 20.-1 = 2. r..(— j (2.20) 

This result looks very nice when translated in terms of turbulence, since it 
predicts (using (2.15)) that the difference of velocity between two points at 
distance r will scale like r'^K/CK-i _ g^j^^ g^jg^ ^j-^^^ time scales and length 
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scales are related through t'^^ ~ r, with C,k = 3/2 (Richardson diffusion). 
Note that r* is simply the convective time across the injection length A. 

These results thus precisely reproduce the Kolmogorov scaling, which is 
here derived from a very simple argument on the directed polymer problem! 
However, it turns out that this result is wrong, because we must go beyond 
the linear approximation. 

Technically the reason is in the integration which leads from the expres- 
sion (2.18) of the average position in Fourier space to the scaling expression 
(2.19). Clearly the corresponding integral over frequencies u is divergent at 
small uj, and it turns out that the result is proportional to (t — t')^T,where 
T is the total length of the polymer. One sees that, even if t — t' < t*, the 
scaling depends on time scales which are larger than A. Therefore one cannot 
work out the scaling behaviour at small length and time scales within the 
linear - random force - approximation. In the polymer problem one must 
study the case of a random potential, which is a non linear problem with 



many metastable states pT], Translated into the turbulence language, 
these correspond to intermittency effects which are crucial and cannot be 
neglected. The full solution derived in the next section indeed finds a Kol- 
mogorov scaling, but only for the third moment of the velocity difference 
{v{x,t) — v{0,t)Y — The other moments do not agree since all the 
moments {v{x, t) — v{0, t)y with q larger or equal to 1 scale like \x\. Further- 
more, length and time are related by a convective scaling, corresponding to 
the dynamical exponent z = ^ = 1. 



3 The replica variational approach 
3.1 The replica solution 



We now turn to the problem of computing the average 'free energy' log Z(x, t) 
of the directed polymer over the random forcing V{x,t). The procedure we 
use is standard, and we follow closely ^]: we first express log Z{x,t) 
as the 'zero replica' limit logZ = lim„^o , average which generates 
an effective attraction between replicas, which we treat using a Gaussian 
variational Ansatz. The quality of this approximation, and the regimes of 
dimension and Reynolds number in which it becomes exact, will be discussed 
in section (3.3). 

We proceed as usual by writing the average of as the partition function 
of the n— replica Hamiltonian Tin, which reads: 



TCr 



1 " /■* 
-T dr 



dr 



WN 



Ej^dre^p 

a,b 



[XaiT 



2iVA2 



(3.1) 
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where we have added a 'mass' term for regularizing intermediate computa- 
tions, which we will set eventually to zero. We furthermore consider periodic 
polymers, for which x{t) = x{0). This is not exactly the same problem as the 
one with free ends which is in direct correspondence with Burgers's equation. 
We shall first work out the periodic case, and in section (3.3) we discuss the 
changes for the free case. 

In order to handle the problems of the metastable states, the idea is 
to use a variational method and approximate 7i„ by an effective Gaussian 
Hamiltonian Tiy, which we write in Fourier space as: 

^v^l^Jdu^ Xa{-^)G-Jiu;)x,{^) (3.2) 

a,b 

where / du stands for Note that is isotropic in real space; however 

its structure in 'replica' space is arbitrary. The trial free-energy obtained with 
depends on Gab and reads: 

N 

=< Tin >v - yTr log G (3.3) 

where < ... >^ means averaging with the Bolzmann weight associated to the 
trial Hamiltonian Ti,^. Remember that the temperature scale has been set to 
one. The calculation of < Hn >v using Gaussian integrals is straightforward 
and leads to: 

<nn>.= ^(Y.jMcUj' + f^)Gaa{u;)-Wj:(l + ^) M (3.4) 
\ ^ oh J 

where T is the total length of the polymer and Bab = J duj[Gaa{^) +Gbb{uj) — 
2Gab{i^)]- The variational statement is that is greater than or equal to 

the 'true' free-energy. We thus look at the optimal Gab such that ^^"I'^l ^ 0. 
This leads to the following self-consistent equations: 

.-11 ^ ^ W f, Bab' 



and 



[G-'Uiu;) = -^[^ + J^,) (a^b) (3.5) 



Gaaiu;) + Gabioo) = G,{u) = — (3.6) 

The task is now to solve these equations using some Ansatz on the structure 
of Gab- This has been discussed in full details in ||2l]]. We keep here to 



dimensions N > 2. There are two regimes of Reynolds number, separated by 
a critical value Rcc = [2(1 - 2/A^)(i-^/2)]i/3_ p^^ ^ ^^^^ ^^le solution is a 
'rephca symmetric' one with Gab{^^) = 6abGc{uj)+G{u!). The two propagators 
Gc and G are easily computed. Translated in terms of the velocity, this just 
corresponds, for an infinite size system (but for a finite injection length A) 
to a vanishing velocity field. For a finite 'box' of length L = ./^, one finds 
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that V ~ ^^2° , where xq is a random (Gaussian) time dependent variable 
depending on the forcing history. When one increases the forcing beyond 
Rcc the correct solution is the so-called 'one step' replica symmetry breaking 
scheme. This amounts to parametrizing the off-diagonal elements Gabi^^) 
with two functions Go{uj),Gi{lj), depending on whether the replica indices 
belong to the same 'block' or to different blocks. The size of these blocks are 
furthermore parametrized by a number m ranging between and 1, fixed by 
imposing that — " ^ '"^^ = 0. We refer the reader to |2T]| for a more 
detailed discussion of this construction, but will recall below its physical 
interpretation, on which we shall heavily rely to discuss our results in the 
turbulence language. The relevant formulae for inverting such matrices are 
given, for completeness, in Appendix A, together with the basic calculation 
steps needed to solve Eqs. (3.5-6). 

One important final equation is the one fixing m. We find that: 

(1 - j^) ' = 2m^—— = -m^Re^ (3.7) 

where we have used the definition of the Reynolds number (Eq. (2.4)). Note 
that Eq. (3.7) is well behaved in the limit A^ oo, where it becomes 
2gm — [jjiReY. We shall see later that when Re > Rcc, the breaking of 
replica symmetry corresponds to a non trivial structuration of the flow into 
large size structures which are'cells' of size A, inside which the velocity is of 
order va ~ R^'^ (independently of the box size L). We shall thus identify 
Rcc as the critical Reynolds number for the onset of turbulence. 



3.2 Physical description of the solution 

Let us see how this is encoded in the one-step solution, by first looking at 
equal time correlations. Within the Gaussian variational Ansatz, one finds 
that 

r[{Ut)}] = ^"'Eexp i-lT.lQ'%-)Mh)Ut)xb{t)] (3.8) 

where Qab = J duG abif^j) -i and X^tt denotes the sum over all the n\ permuta- 
tions of replica indices. Taking the limit where the mass term is zero, we 
find that (see Appendix A): 

Qo ~ /x^/'-' ^ 0, Qi-Qo^oo, Q-Qi= ,!^^ln. (3-9) 

U N ) 

In the very large Reynolds limit, on which we focus now, the parameter 
m scales as 1/Re and one finds Q — Qi — mA^. 

The physical interpretation of the replica probability distribution Eq. 



(3.8) has been worked out in and is particularly simple in the present case 
where Qq is zero and Qi~Qo is very large. Remember that one is dealing with 
a polymer in a random environment. For each 'sample' VL (i.e for a particular 
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realisation of the disorder V{x,t)) the probabihty distribution Pn{x) for the 
'end point' of the polymer x has a certain shape. The 1 step replica symmetry 
breaking variational approach assumes that this distribution can be written 
as a weighted sum of Gaussians: 

Pnix)=Y.W^exp(^-^[x-f^]^^ (3.10) 

where 6 = Q — Qi, are random weights chosen with a probability propor- 
tional to W~^~"^{1 — W)"^~^ and Tq, are uniformly distributed in the 'box' of 
size L (A more formal description of the solution is presented in Appendix 
B, together with useful technical details). The Wa and the Tq, encode the 
particular features of the sample fi, and must thus eventually be averaged 
over. 

The above construction was restricted to a certain 'time' (or length of 
the polymer). Similar considerations also enable us to construct two-time 
correlation functions, and the result is very simple: it amounts to let the 
acquire a time dependence. More precisely, Tq, are independent Gaussian, 
time dependent variables such that (again taking the /i = limit): 

C{t - t') ^ < [f^it) - Ut'W > = - I d J^-'^'f-''^^ , (3.11) 

where S = l/(4c(Q — Qi)^)- Eq- (3.11) shows that a characteristic time 
scale appears, given by = 2^^^-^, which is nothing but the 'Larkin- 
Ovchinnikov' time r* encountered above and corresponds to the convective 
time across the injection length. For r ^ r*, one finds that C(r) = 23A'^^, 
whereas for r ^ r*, one finds C(r) = |A^(;^)^. 

Eqs. (3.10-11), together with the value of m and 5, are the central re- 
sults of this paper, from which we shall derive in the next section the very 
interesting statistical properties of the velocity field, in particular the exact 
calculation of the full probability distribution of v{x,t) — v{x',t). Before 
turning to this computation, we first discuss the validity of this replica solu- 
tion. 



3.3 Discussion of the replica solution 

One is not able to solve exactly the directed polymer problem. The vari- 
ational method to which we have resorted is however known to be a good 
approximation to the real behaviour of the directed polymer (|^, |3^). An 



interesting advantage of this approach is that it gives exact answers for the 
thermodjTiamics in the limit of a large number of dimensions N —>■ oo. Tech- 



nically, as shown in , this comes from the fact that this variational Ansatz 
can be seen as a resummation of the 'Hartree' diagrams in the perturbation 
theory for the correlation function, which are the only ones surviving in the 
large N limit. However beyond the thermodynamic potentials, there are 
some quantities for which this approach fails to give the right answer, even 
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for N ^ oo. It is unfortunately the case of the velocity correlation func- 
tion defined in (2.15), for which our calculation may only become exact in 
the large N and large Reynolds limit. The problem can be understood as 
follows: the replica calculation is expected to reproduce faithfully the ther- 
modynamical behaviour of the system, by correctly describing the low-lying 
energy states. The weight of these low-lying states is found to be distributed 
as oc Vr~^"™(l — W)"^~^. Hence the moments [Pn(x)]'' which are primarily 
determined by the low-lying states (i.e. by those with large weights) are 
expected to be accurate, while those which are sensitive to the W ^ part 
of the distribution are not. This is the case when 77 < m, in particular for 
logPn{x) which we need in order to calculate the velocity in the Burgers 
language. Hence, only when m — 0, i.e. at infinite Reynolds number, can 
the method be reliably used to determine exactly the statistical properties of 
the velocity field. To summarize this discussion we expect that this approach 
may be exact only when N ^ oc and m —>■ 0. This issue can be studied 
from the solution by a study of the correlation identities, some of which will 
be checked in section (4.3). In any case the picture which emerges from the 
discussion in the next section is an appealing one which seems to be a good 
approximation even at small A^. 

Eventually, we want to discuss an important detail. The previous replica 
calculation was performed, for simplicity, in the case of a periodic polymer, 
i.e. one for which x(t) = x{0). From the Burgers equation point of view, 
however, the end point of the 'polymer' must be left free. As argued by 
Nelson and Vinokur in a different context , the difference between points 



far 'inside' the chain and boundary points is the fact that the statistical 
weight of the former is given, in quantum mechanical language |]3^, by the 



product of two propagators, while the end point only needs one. If the ground 
state wave function of the corresponding quantum mechanical Hamiltonian 
decays exponentially on a characteristic length scale i, this means that the 
lateral fluctuations of the polymer's end point decay on length scale i, while 
the fluctuations for points in the bulk decay on a length scale £/2. The 
study of the system with periodic boundary conditions which we performed 
in Sect. (3.1) only deals with these fluctuations in the bulk. The simple 
consequence of this analysis for our problem is that formula (3.10) can be 
used for free end points as well, although the value of S must be doubled, 
i.e.: 

S = 2mA2 . (3.12) 



4 Structure of the velocity field 
4.1 Qualitative Arguments 

We can now use our 'dictionary' between the two problems, in particular the 
fact that the velocity field is given by the derivative of the free-energy of the 
polymer problem with respect to x. For a given realisation of the forcing. 
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this free-energy is simply the logarithm of the probability distribution Pn{x) 
given in Eq. (3.10); one can thus express v as: 



We shall study the statistical properties of this velocity field at large Reynolds 
numbers. We recall that from the previous analysis, we have S = 2mA^ ~ 
2^^^A'^/Re. Typical snapshots of such a velocity field in A'^ = 1 or = 2 
dimensions are given in Fig. 1 a,b. Very clearly, a cellular structure appears. 
Within each cell, the velocity field is radial v ~ ^{x — fa), with a rapid 
variation (shocks) across the boundaries of these cells. 

This structure can be understood qualitatively from Eq. (4.1), using 
the important property that the weights Wa have a very broad distribution. 
More precisely, we know that Wa = exp(— /q,) / E^ exp(— /^), where the fa are 
independant random variables with an exponential distribution increasing as 
exp(m/). For m small, there is a very strong hierarchy between the smallest 
/ appearing in the sum. For instance the gap between the smallest fa and 
the next scales as 1/m oc i?e. So the low lying fa, which are the only terms 
contributing to the sum, behave at large Re as fa — faR^, and the velocity 
field in this regime takes the form: 



2u 
T 



Ea(^-^a) expi 


'Re 


J a 2^/3 A2 


) 


Ea exp (i?e 


f [x—ra)'^ 
J a 2V3A2 


) 



(4.2) 



For a given x and a large Re, the major contribution to the sum comes 
from the largest term. This dominant term will suddenly switch from -say- 
a to f3 when the space dependent Gaussian factor compensates the weight 
difference, i.e. when {x* — faY ~ (^ ~ fpY + m^ This leads to a typical 
size of cells of order oc A, i.e. the injection length itself. The width £ 
of the shock separating two cells is obtained by writing that ~ i-^- 

£ oc -C A, which shows that the notion of cells is indeed well defined at 
large Reynolds. This construction allows one to guess the structure of the 
probability distribution of the velocity difference between nearby points. Let 
us present the argument in = 1 dimension (its generalisation to higher 
dimensions is straightforward), and call x and x + r these two points. The 
probability that a cell wall is present within the interval [x, x + r] is obviously 
given by p = In the limit where p is small, the velocity difference u — v{x-\- 
r) — v{x) is equal to ^ with probability 1 — p, and of order ~ ^ = i>a 
with probability p. Hence: 

P,M~(l-rt.(»-H^)+A/(ii) (4.3) 

where /(.) is a certain scaling function. This is true in the limit Re ^ oo 
where the internal structure of the shocks can be neglected. We will see 



12 



below that precisely such a form is obtained from an exact calculation. Let 
us now estimate the various moments of u using expression (4.2): 

'^-il-p)C-£y+pA,vl, (4.4) 

where Aq = J duu'^f{u). In the limit where £ ^ r ^ A, one finds that for 
all q larger than 1, one has 

^ = Aqvl^, (4.5) 

whereas for q < 1, one finds 

T^^vii^r. (4.6) 

The velocity field is thus strongly intermittent, with a 'multifractal' spec- 
trum given in Fig. 2. Note however that the q = 3 moment scales a la 
Kolmogorov, i.e as in the real turbulence problem. The presence of large 
scale structures (shocks) forming an N— dimensional froth-like pattern (see 
Fig. 1-b) is responsible for such a strong intermittency. Large scale structures 
in 'true' turbulence are similarly thought to be the origin of the experimen- 
tally observed intermittency, which is however much milder (see Fig. 2). The 
deep reason of this difference is probably related to the 'dimension' of the 
large scale singularities, which is iV — 1 in the present case and only 1 for 
vortex lines in hydrodynamical turbulence. 

This strong intermittency, due to the cell and shock structure, has already 
been discussed in the one dimensional Burgers turbulence in the decaying 
(unforced) case (||3^, ^). To the best of our knowledge the forced case 
has not been studied. Furthermore we are able to derive the exact form of 
the probabihty distribution function of velocity differences, to which we now 
proceed. 

4.2 Exact results at infinite Reynolds number 

Using the direct evaluation of the moments of the velocity difference from Eq. 
(4.1), one can establish, after rather long manipulations detailed in Appendix 
B, the following form for the full distribution of longitudinal velocity differ- 
ences. Let us first start by the equal time case, and write u = v^{x) — v^{y), 
where the component '1' is along the x — y axis, and / = (The other 

components would be treated in a similar manner). The result then reads: 




Pi{u)=6\u-l—^-j e^' j_^DhR{h,\l\) 

/oo roo POO I Oy I \ 

dh Dt Ds6 \u+ {t + s) — J— sign{l)\ R{h, 

-oo J-h-\l\/2 Jh-\l\/2 \ ^ m / 



(4.7) 



13 



where R{h, I) = le-'^^/'^Moih - 1/2) + e^^/'^Moi-h - \/2)] ^ {Mq is related 



,2 



to the standard error function - see Appendix B), and Dx = -j==e 2 . Note 



2-K 

that Eq. (4.6) has precisely the structure which we guessed using qualita- 
tive arguments (Eq. (4.3)). We have tested this formula numerically by 
generating a one- dimensional velocity field using eq. (4.1), for Re = 100, 
and directly computed the second moment as a function of p, which we 
compare in Fig 3 with the exact formula for Re 00 obtained by inte- 
grating (4.7). For <^ p ^ A, corresponding to the 'inertial range', one 
finds = ^^^^^ ^^^^^ , or in terms of the 'traditional' energy spectrum 

E{k) = k^-^{v{k)v{-k)), E{k) oc k'"^ for all A^. Note the saturation for 
distances much larger than the cell size A, on which we shall comment later. 
We have also computed the weight of the 5 peak and compared it to the one 
obtained analytically, with good agreement again. 

We have also extended our analysis to two interesting situations. First of 
all, one may consider the case of a finite Reynolds number. The calculation 
of the second moment v? shows that a new length scale appears^ = 
separating the above linear regime of v? at larger distances, from the regular, 
quadratic behaviour at small length scales: 

^ = — r^; r < ^f (4.8) 
2v ^ ' 

which is the standard result (up to the numerical prefactor) obtained in tur- 
bulence for length scales smaller than the dissipation length (see e.g. 



Note that the 'dissipation length' Id depends on the moment of the velocity 
one wishes to calculate: this is another consequence of intermittency. Match- 
ing the regular behaviour (4.7) with (4.4) or (4.6) suggests that the (^^ order 
dissipation length scales as = ARe^-^i ; only for g = 3 does one recover 
the usual Kolmogorov dissipation length i^^^ = -j^ji- 

From our Ansatz for the velocity field (4.1), and the fact that fa{t) evolve 
according to (3.11), one can see that the full velocity correlation function has 
the following scaling form (in the limit ^ r ^ A): 



[v{x + f,t)-v{x,t'W = r g\ . (4.9) 

l2C{t'-t) 

where C(r) is defined in Eq. (3.11) and g{.) is a certain scaling function. 
We have obtained the precise form of g in the case where t and t' are well 
'inside' the polymer, but did not attempt to compute it for 'end points' - 
where additional numerical factors would appear (see the discussion at the 
end of section 3.3). 



In the limit t-t' < r* (see Eq. {?>.ll)), ^2C{t - t') = v/^\t-t'\: Eq. (4.9) 
then means that fluctuations travel in a ballistic way with a 'velocity' fixed by 



■^Note that I'^/p is the width of the shocks t introduced above 
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the injection scale. For larger time scales, the dynamics recovers a diffusive 



character, since ^2C{t — t') oc Ay^^— ^. Note however that the effective 

diffusion constant i^turft = ^ oc i?e z/ is enormously enhanced compared to 
its 'bare' value - as is the case for usual turbulence. Prom Eq. (4.9) in 
the limit r — 0, one finds that for coinciding points, the velocity difference 
grows with time as: 

[v(f , t) - v{x, t')f oc ^2C{t - t') (4.10) 

4.3 'Sum-Rules' and the large Reynolds limit 

There are a number of 'sum-rules' that the correct solution of the Burgers 
equation should satisfy, allowing us to test our prediction for the structure 
of the velocity field, as given by Eq. (4.1). The first and most interesting 
one physically is the 'energy conservation', i.e. 



IdiP 



2 dt 



0^ -V- {v-V)v + u V- Av + v- f (4.11) 



The last term of this expression is the injected energy, and is equal to +^6. 
The two first terms represent the 'dissipated' energy - notice that contrarily 
to the incompressible Navier-Stokes case, this dissipation is non local (i.e. 
it cannot be expressed in terms of velocity derivatives only). We find (see 
Appendix C, Eq. C.ll), that the total dissipation is given by — ^€(1 — m). 
Hence we find that the energy conservation is indeed exactly satisfied in the 
limit m — > {Re — > oo) at fixed e (i.e. for A — > oo or ^ 0). Interestingly, 
in this limit, this sum rule is true independently of N. Eq. (4.11) can also be 
interpreted slightly differently. Suppose that at t = the forcing is switched 
off. The subsequent evolution of the energy density is then given by the 
dissipative terms, i.e. iVA^/^^^^ ~ leading to e ~ or else for the 

velocity scale v^it) ~ y. 

We have also checked that ^^^^ = in the hmit it!e — > oo for all N, while 



— identically for all and Re. The calculation of velocity correlations 

dt 



for non coinciding points, such as ^v{x)v{y), are much more intricate, and 
are currently under investigation. 

Hence these direct checks of certain (local) correlation identities suggest 
that our Ansatz might be exact for large Re. 



5 Discussion 

There are quite a number of points worth discussing under the light of the 
previous results, before commenting on the similarities and differences with 
the 'real' turbulence problem. 
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5.1 Burgers turbulence in small dimensions 



As stated above, our results are a priori only exact for large dimensions and 
Reynolds number, the method we have used becoming a variational approx- 
imation in finite dimensions. By comparison with other problems where this 
approach was used, however, one expects that our solution describes faithfully 
the physical situation, even for relatively small values of N. Furthermore, all 
the sum-rules which we have checked so far (Section 4.3) are satisfied for any 
N in the limit Re — oo. Hence, we believe that the forced Burgers equation 
in say = 3 dimensions will produce a cellular arrangement of the fiow 
pattern, as described above. It would be extremely interesting to test this 
prediction numerically, as well as to measure the intermittency corrections 
- which might be weaker in finite dimensions than the ones obtained within 
our Ansatz. 

For N <2, one however knows from previous studies that the structure of 
the replica solution changes: instead of a 'one step' breaking, a full continuous 
breaking scheme is needed. This raises the interesting question of knowing 
whether such a scheme could describe a more complicated behaviour of the 
moments of the velocity field than the one encountered above. We do not have 
a complete answer to this question, but it seems that the solution obtained 
for N <2 does not lead to a qualitatively different picture, the reason being 
that all the length scales appearing in recursive construction of the sample 
dependent measure ||21[] (i.e. the width of the Gaussians), as well as the 
'm' parameters describing the weight distributions, are all of the same order 
of magnitude. Thus, forced Burgers turbulence in one dimension is expected 
to look very much like the picture shown above (Fig 1-a), as indeed is well 
known from studies of 'decaying' turbulence from random initial conditions 

More work on this aspect would certainly be interesting, in particular to 
study a variant of the present problem for A^ = 1, where / - rather that - 
is a random noise (see |38|). In this case also, a full breaking of replica 



symetry is needed to describe the velocity field. 



5.2 Long-distance scaling of the velocity field 

As discussed in Section 2, the main difference between this work and previous 
studies on the directed polymer problem lies in the regime of length scales. 
Usually, in the directed polymer problem, one assumes that the correlation 
length of the potential A is very small, and one is interested in the long- 
distance X ^ A scaling behaviour of - say, the free-energy, where non trivial 
exponents appear. On the other hand, in line with most studies of the 
turbulence problem, we focused on the contrary on the small scales x <^ A 
regime, where an inertial range appears, characterized by an energy cascade 
E{k) oc A;~^. The existence of these two regimes for the Burgers equation 
suggests that a similar situation may also occur in 'real' turbulence, where a 
non trivial scaling regime could exist for x ^ A, characterized by a velocity 
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correlation converging towards its asymptotic value as a power law: 





+ f) - 2t;(f)2 - — — ^ + ... r>A (5.1) 



-2 



or equivalently, by a new exponent for the energy spectrum: E{k) oc 
where u and ( are the generalisation of the exponents defined for the poly- 
mer problem in Eqs. (2.13-14), corresponding now to the usual large scale 
situation r ^ A studied in directed polymers. Numerical results on directed 
polymers in 3-1-1 dimensions Jl^ give C ~ 0.6, u; = 2^ — 1 ~ 0.2. It would 
be very interesting to analyze experimental data beyond the injection length 
A along these lines. 



5.3 Passive scalar dispersion 

A subject of recent debate is the behaviour of a passive scalar in turbulent 
flows, which also shows experimentally 'anomalous' density fluctuations - 
in particular exponential tails Although a detailed study is beyond 

the scope of the present paper, it is interesting to discuss qualitatively this 
problem for Burgers' turbulence. Let us suppose that the diffusing scalar 
obeys a Langevin equation of the form 

^ = v{x,t) + fj{t) (5.2) 

where v{x, t) obeys the forced Burgers equation Eq. (2.1), and ri(t) describes 
the molecular diffusion, with diffusion constant equal to D. In the 'adiabatic' 
limit where the evolution of the velocity field can be assumed to be small 
compared to the equilibration time (i.e D ^ ^aA), the tracer concentration 
6{x,t) is simply given by the Boltzman equilibrium: 

^(f,t) ocexp-^^^^ (5.3) 

where h{x, t) is the velocity potential introduced in Eq. (2.5). Using h{x, t) = 
21/ log Pq{x , t) + est., one finds that: 

9{x,t)^[Pnix,t)]~"-^ (5.4) 

Now the full distribution function of Pq can be exactly calculated, and 
is found to be a totally asymmetric Levy distribution with index /i = m. 
This in turn allows one to obtain the full distribution of 6, as 

'p^0)=0~^-^LUO'^) (5.5) 



Hence, we find that the tailoi the tracer density distribution is a power-law for 
small 6, V{6) ~ 6^^^~ , which becomes broader and broader as the Reynolds 
number increases (m oc Re~^). For large 9, V{9) ~ 9^^^ exp 



-cst^^™ 
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Note however that our adiabatic assumption breaks down precisely when 
— < 1. It would be interesting to understand the nature of the fluctuations 
of Q in the large Reynolds regime. 

Hence we have argued that the fluctuations of tracer concentration are 
strongly anomalous in Burgers's turbulence. In more physical terms, Eq. 
(5.2) shows that the tracer is convected towards the shock regions, where 
the concentration piles up. High density 'sheets' of particles spontaneously 
form, revealing the 'froth-like' structure of the flow. This is similar to what 
happens for the Burgers equation with random initial conditions (but no 
forcing) and might be relevant for astrophysical applications 



Finally, it must be noted that the temporal dispersion of tracers is convec- 
tive at short times, since it essentially follows the velocity fleld which evolves 
according to Eq.(3.11). This is, as noted above, in contrast with Richard- 

3 

son's diffusion, which states that x . Convective dispersion is however 
rather commonly observed (see [Q). 

5.4 Turbulence versus 'Burgulence' 

The differences between Burgers' turbulence (coined 'Burgulence' by Fournier 
and Frisch in [Q) and hydrodynamical turbulence have been discussed many 
times. The most important one is the absence of vortex line singularities in 
a potential flow, which are thought to play an important role in turbulence 
^ . In our case, singularities are concentrated on — 1 structures, giving 
rise to much larger 'intermittency' corrections. 

Second, the energy is only dissipated by viscosity at small length scales 
in turbulence, while there is an additional dissipation term in Burgulence 
(which is, as discussed above, non-local). This might be another important 
difference although, as shown above, it does not prevent the existence of a 
well-deflned inertial range where the energy spectrum follows a decay. 

A third difference invoked by Kraichnan is that in fluids, the incompress- 
ibility condition is maintained by the underlying pressure fleld, which would 
play an important role in the dynamics of the fluctuations and be at the ori- 
gin of the Richardson-Kolmogorov scaling x ~ t^^^, instead of the convection 
law X ~ VAt found here. On the other hand, the Burgers and incompressible 
Navier-Stokes equations look very similar, in particular from a dimensional 
analysis point of view, which is at the heart of Kolmogorov's argument; fur- 
thermore, the third moment of the velocity difference indeed scale in the 
same way (at least for large N). 

On a technical level, the simplest closure scheme for turbulence is Kraich- 
nan's DIA, which leads both for the Navier-Stokes equation and for the Burg- 
ers' equation, to a convective dynamics x ~ v^t and to a k~^^'^ energy spec- 
trum. However, DIA's extension to Lagrangian coordinates (LDIA) was ar- 
gued by Kraichnan to reproduce exactly Kolmogorov's scaling for the 
Navier-Stokes equation, but not for the Burgers case, where results similar 
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to some of ours are obtained Q It would be interesting to investigate the 
precise relation between the two seemingly very different approaches, and in 
particular to compare LDIA to our results. 

5.5 Conclusions 

We have thus used a method inspired from spin-glasses to investigate a toy 
model of turbulence, and to propose an Ansatz for the structure of the ve- 
locity field, which should become exact in high dimensions and when the 
Reynolds number is also large. We find that beyond a critical Reynolds 
number, there exists a well defined inertial range where the energy spectrum 
decreases as The scaling variable is where va is the velocity at 

the injection scale A. The third moment of the velocity difference scales 
linearly with distance, a la Kolmogorov, but strong intermittent corrections 
come into play, due to the presence of shocks localized on a froth-like, cellular 
pattern. The full distribution for the velocity difference is obtained exactly. 
Interesting scaling results are also argued to hold at scales larger than the in- 
jection scale, and we suggest that experimental data on grid turbulence could 
be analyzed accordingly. Our results are presumably qualitatively correct in 
low dimension; numerical simulations would be welcome. We have discussed 
qualitatively the passive scalar problem and we have found power-law tails in 
the concentration distribution, reflecting the localisation of the tracers near 
the shocks. 

From a technical point of view, it would be interesting to understand the 
precise relation between the present approach and Kraichnan's Lagrangian 
DIA, which gives the same scaling as the ones obtained here. From a different 
point of view, one could also generalize this work to the case of a space- 
correlated forcing term f{x, t) with a power-law correlation function. This is 
the starting point of the renormalisation group (RG) analysis of turbulence: 
the exponent describing the decay of the power-law function is chosen as to 
reproduce Kolmogorov's scaling and an RG procedure |^ is applied to obtain 
adimensional prefactors (the Kolmogorov constant). The model studied here 
could provide an interesting benchmark to discuss the validity of such a 
procedure. In fact it can be used to test all the various approximate methods 
which have been introduced in the study of fully developped turbulence. 

Finally, our Ansatz for the velocity field [Eq(4.1)], which is inspired from 
our 'replica' approach to the problem, has some interesting mathematical 
properties (see in particular Appendices B and C). It could be fruitful to 
generalize Eq. (4.1) to describe rotational flows. 
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Appendix A: Solution of the saddle-point equations. 



The starting point of the one-step replica symmetry breaking calculation 
is the expression of the free-energy in terms of G{u}),Gq{uj),Gi{uj) and m. 
Using the expression (AIL 11) given in |2l[] for the trace of the logarithm of 
a one-step Parisi matrix, we find that: 



'N/2 



+ m(l + 



Bo . 



-N/2 



} 



-- / du. 



m 



■log(G'(cj) -mCoH - (1 



Ga(uj) 



G{uj) - mGoiuj) - (1 - m)Gi{uj) 



1 — m 



m 



log(GH-GiH) 



(Al) 



where Bi = 2J dLj[G{uj) — (We have set the temperature to T = 1). 

Differentating with respect to G^u), Go{uj), Gi{uj) yields equations (3.5-6) 
specialized to the one-step solution, with 



Gciuj) = G{uj) - mGoico) - (1 - m)Gi(cj) 



G,(u) 



and 



G^iu) - Goiuj) 



1 S 



(A2), 



(A.3) 



(A.4) 



(/i + ccu^) m(/i + cuj'^ + S) 

where we have introduced S = m{[G^^]o — [G^^]i). Using the definitions of 
Bq, Bi, we obtain the equations: 



Bi-Bq = duj 



m 



S 



H + au"^ yU + cu"^ + S 



{A.5) 



and, using (A. 2) 



mBo+(l-m)Bi = 2 / duj[G(uj) -mGoiuj) - (1 -m)Gi(uj)] = 2 / du 

J J jl + OjJ^ 

{A.6) 

from which we can deduce Bq and Bi as a function of S. Inserting these 
values into Eq. (3.5) then leads, in the limit fi ^ 0, to: 



[G-%r^O; S=~m[G-']i 



mW 



[1 + 



1 



-l-N/2 



(A7) 



Differentiating now (A.l) with respect to m, after a few manipulations, leads 
to 



- = 2m^W\l + 



-N/2 
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Eqs (A. 7-8) allows us to obtain both S and m as: 
and 

S = ^ — (A.IO) 
Using Eqs (A2-7), we thus obtain in the limit // — > 0: 

Qc=7r^, Qi-Qo^TT^^, g-gi = -i= (All) 
and Qo //^~ ■ Prom this, we also obtain the inverse R of the matrix Q as: 



D ^ n I D ^0 4 
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Appendix B: Probability distribution of the velocity difference 



In this Appendix we compute the probability distribution function (pdf ) 
of the difference of velocities between two points x and y, projected onto 
X — y, at large Reynolds numbers. We shall proceed in three steps: first we 
show how to compute the first few moments for general Reynolds numbers. 
Then we show how the expression simplifies for large Re. This simplification 
is such that we can extract in this large Re limit all moments and deduce 
from it the pdf. 

Our starting point is the result from the replica computation, which pro- 
vides the following random process to build up the velocity distribution: at 
a given time, the velocity field is given by: 

.... 2i^ EaW^a(:g-r-;)e-("^-W2^ 

The sum over a goes from 1 to M. The points are uniformly distributed 
in the volume [—L/2, The weights Wa are random numbers drawn as 

follows: First one chooses M "energies" fa at random, they are identically 
distributed independent random variables with a probability distribution: 

V{{) = t^^^^-^M{i-{) {B.2) 
Then the weights are given by: 

We shall let M, L, and fc go to infinity together, keeping the density of 
states Me~"^^'' / L fixed. The and are uncorrelated. We shall denote 
respectively by Er and Ef the expectation values with respect to these two 
sets of random variables. 

Let us first evaluate the second moment. We take x — y along the first 
axis, and compute the longitudinal correlation: 

vi(x)vi(y) = (^y) 




e-fc-f^e-[(^-''-)^+(y-^/3)^l/25 



[BA) 



We shall compute separately the two contributions to this correlation which 
come from the terms a = (3 and a ^ (3 in the numerator of (B.4). We thus 
write: 



vi{x)vi{y) = l^Y J (911 +912): (B.b) 
where gu and gu respectively contain the terms a — (3 and a ^ (3 in (B.4). 
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To compute gn we write the denominator of (B.4) in integral form, which 
gives: 



poo roo 

gn = MEr,f / dA / dfxe-''^"{yi - r«)i(y - f„)ie- 
Jo Jo 



[(x-r,)2+(y-ra)2l/26 



exp 



n (exp 



As for the mean over the free energies, we shall use repeatedly in this Ap- 
pendix the following formula, valid for large fc- 



After taking this average over the free energie, (B.6) becomes: 
Jo Jo 



n [l - e-'"^=r(l - m; 



(S.8) 



We now average over the values of fa and fp, with a uniform measure in 
a box of size L^. In the limit of large L, M, fc at fixed density we get: 

Me~^^'^ /■oo /■oo 

gn — mT(2 — m) / d\ j dji 

L Jo Jo 



j df{x — f)i{y 



2/oA\™-2 



exp 



(5.9) 



It is convenient to rewrite — > A//', and integrate over A, which gives: 

roo 

gu = (1 - m) / d/jL 
Jo 



(5.10) 

(In this expression the / dr is over a single r variable, namely the compo- 
nent of r in the direction of x — y, the other components have already been 
integrated out, their contributions cancel between the numerator and the 
denominator. We denote x — Xi, and similarly y — yi) 
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It is straightforward to perform the same steps for the second contribution 
gi2 in (B.5). One finds: 



m—l 



J 



m—l 



j dr (1 



g (x-r)V25 ^ ^g-(2/-r)V25y 



-2 



(5.11) 



It should be clear that these techniques allow to write the velocity correlations 
of low order in a relatively closed form (this means reduced to some finite 
dimensional integrals, where all the averages over r and / have been taken 
care of, as well as the limits M — > 00, L — > 00 and fc — > 00). However 
the expressions are comphcated enough, especially when one goes to high 
moments, and we have not found a closed form for the pdf in general. 

Fortunately the situation simplifies in the limit of large Reynolds num- 
bers. Let us first work out the expression of the second moment in this 
Re 00 limit. We know from the replica solution that in this limit the 
width 5 scales as 5 ~ 2m A^, where A is the scale at which energy is injected, 
and m, the breakpoint in Parisi's RSB solution, behaves as m ~ 2^^^ /Re. 
Hereafter we shall use m instead of Re, and we want to understand the small 
m limit of the velocity pdf. 

We proceed and first work out the small m limit of the two point cor- 
relation. In gii we change variable to x — XyJS/m, y — y^S/m; We shall 
take the m ^ (large Re) limit keeping x, and y fixed. The algebraic dis- 
tance between the two points is measured hy I = (x — y) = {x — y)i/ (Aa/2). 
We also change the dummy integration variables in (B.IO) from r to ^ = 

r^J 5/m — {x + y)/2, and from ^ to h = — mln(/i)/|/|. This gives: 



9ii 



[l-rn] 



dhe-^\^\''^ 



where we have introduced the functions (defined for integer ki, k2): 

dz 



(5.12) 



KkM^ ljL{l/2-z)'^{-l/2-z) 



k2 



27r 

-z'^/2+lz{ki-k2)/2m ^^lz/2m 



+ e 



-h\l\/m—lz/2m 



m—{ki+k2) 



(5.13) 



Their small m limit is easily worked out by a saddle point integration. The 
result for /ci > 1, A;2 > 1 is: 



Kk2ih) -m^o ^sign(/)'=^+'=^ e'^\im2/m-i/2) e 



-h^/2 



27r 



k2 



+ h 



r(A;i)r(A:2) 
r(A;i + k2) 



[B.U) 
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While the results for /ci/c2 = read: 



(i:i(iis/$-oo) 



/o,o(/i) 



.0 e 



fi/S-h\l\/2 



In these equations we have used the definition: 



(S.15) 
(E.16) 



Therefore we get for the m — > limit of gu: 
gii^ I dhe-^^'"^ 



-f-f 

m\ TT Jo 



{B.17) 



'e-'^mM>{{ - f ) + y\i\/^M,{-{ - f )" 

The small m limit of the second contribution to this second moment can be 
worked out with the same technique. One gets: 

^ r 

912^ / dhe 

m TT Jo 



1 



=^2/4 



(B.IS) 



We remember that from (B.5) the velocity two point correlations between 
points at distance lA\/2 equals (2iy/5)'^{gu + (722) • It is therefore finite 
when m ^ 0. It is interesting to work out its small / limit, which gives the 
correlation between two points whose distance is small with respect to the 
injection scale A but large with respect to the dissipation scale (since we 
have taken the limit it!e — > 00 -or m — > 0- first): 



< vi{x)vi{y) >: 



x-y 



{B.19) 



We now turn to the computation of higher order moments. We are inter- 
ested in arbitrary moments of the velocity difference v{x) — v{y) projected 
onto the direction oi x — y (chosen as the '1' direction. The starting point is 
analogous to (B.4): 



Vi{x)Pvi{y)p' = ( 



YP+P 



ai,...,ap /3i,...,/3p 



^-(f„,+...+f„p+f5,+...+f;3p,) 
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(f - r;ji...(f - r; Jiexp 



{X-Ta^f + ••• + {x-faf 



-fc-{x-rcf/2S 



13 



[B.20) 



We shall proceed as for the second moment. The first step is to exponentiate 
the denominator in (B.20) using 



ApBp' 



dX 



dfi 



Tip) T{p') 



(5.21) 



In order to perform the averages over the values of the free energies 
and the positions in (B.20) we need to distinguish how many indices in 
ai, ...,ap,Pi, ...,Pp' are distinct one from another. Let us suppose that there 
appear in this sequence k different indices, which we call 71, ...,7^. We shall 
call Qj the number of indices in ai,...,ap which are equal to 7^, and q'j the 
number of indices in j3i, Pp' which are equal to 7^. Such a configuration of 
indices, characterized by the number k G {1,2, ...,p + p'}, and the sequences 
qi, Qk, q[, ...(ik (such that gi + ... + = p, + ... + = and for every j: 
Qj + Qj > 1) appears a certain number of times, which we call C{|'^ ^, , 

in the sum over ai, ...,ap,Pi, ...,Pp of (B.20). For each configuration of in- 
dices, we just do the same transformations as for the second moment. We 
shall not repeat them here, and just give the final result which generalizes 
(B.12): 



vi{x)Pvi{y)p' 



^6' ^m' mJ-00 T{p)r{p 



p'h\l\jm P+p' 



7T E m 



m 



k-l 



k=l 



qi, 



E 



c 



(II) 



,U,,,U' 



■Ik 



n 

1=00 



-(ni + nj-$)Jn„u|(0 



MO 
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Let us now work out the leading behaviour of this expression in the small m 
(large Re) limit. From the integrals Iq.^qrih) we get a factor exp(/i|/| I'j/f^) 



which exactly compensates the explicit e 



-p'h\l\/m 



Using 6 = 2mA^ 



we 



find that the contribution to (B.22) from a given value of k scales as 
{u/5)^^^ [5 / m)~^P^P ^/'^m'^^'^\ A close look at the behaviours of the integrals 
Iq^qi shows that a{k) is zero for k = 1,2, and it is strictly positive for k > 3. 



So we can neglect the terms with k > 3 in (|5.5|) . The term = 1 is eas- 
ily worked out. The only allowed configuration of indices has qi = p and 



p', and its degeneracy is C 



00. For the k = 2 term we keep only 



the leading terms which correspond to gi = 0, g2 = p, q'l = p', g2 = or 
to gi = p, g2 = 0, q[ = 0, g2 = p', which also have a degeneracy 1 (*** a 
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verifier Altogether the leading contribution to this moment at large Re 
is: 



'^Np+p'/_^x{p+p')/2 

' 5 ' ^m' 



sign(/)*'+P' j dh 



+ h 



P' g-(/i2/2+/2/8) 

+ h\ y= R{h, 



27r 
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It is a simple exercise to sum these moments and get the q*'* moment of 
the longitudinal velocity difference: 



dm J 



g-(ftV2+«V8) 



27r 



/oo 
-h- 



dt 



\i\/2 V27r 



r 



|i|/2 V27r 



e-' '''{t + syR{h,\l\f 



where: 



R{h, I) = [e-'^'/'A^o(/i - V2) + e^'^^^Moi-h - 1/2) 



, (i?.24) 



(5.25) 



Under this form the moments can be inverted and wc obtain the explicit 
form of the pdf P{u) of the velocity difference u = ^1(2;) —vi{y) between two 
points at distance x — y = /A\/2 at large Re: 



+ 1^ 



/CO rOO 
dh / 
-00 J—h- 



dt 



r ds 



\i\/2 V27r 



i 

Jh- 



\i\/2 V27r 



S[u+{t + s)^sl^signil)^ R{h,\l\f 



{B.26) 
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Appendix C: Correlation identities 



In this Appendix we check that the solution we have found does sat- 
isfy some necessary identities of correlation functions, in the limit of large 
Reynolds numbers. We shall perform the explicit check in the case of the 
"energy" balance. Prom Burgers' equation one finds: 

1 d- 



A= — vi^di^{vPvP) B = v vP-dP-dPyP C = fpyp , (C.l) 

with a convention of summation on the repeated vector indices jj,, p from 1 
to N which will be used in this whole Appendix. The overbar denotes the 
expectation value with respect to various realizations of the force /. We shall 
compute successively the three contributions A.B (from our solution), and 
C (from a direct computation of the energy injected). In the end we shall 
check that their sum vanishes at large Reynolds number, as implied by the 
stationnarity of the forced flow. 

We use the Ansatz for the velocity given in (4.1) and developpcd in Ap- 
pendix B. It will be useful to introduce somewhat more compact notations 
and define: 

where is a random variable which takes value with a probability 
WQ,exp(— (x — fa)/2S). We shall denote as before by brackets the expec- 
tation values with respect to this process with fixed Wa and fa, while the 
averaging over Wa and fa (corresponding to the overline in (C.l)) will be 
denoted as in Appendix B by Erj. (In terms of the directed polymer, < O > 
denotes a thermal average, and the overline denotes an average over quenched 
disorder). With these definitions we have: 



2u 
1 

dPy'' = y (^SP" - ^(< > - <uP ><u>' . (as) 

which leads to the following expression of A: 



2iy\^ „ r „ „ 1 



A = ( - ) ErJ 



-< >< > + -< u'' >< X u^u^ > 
b 



-< u'' >< X u'' X > 



1 

The same steps give the expression for B: 



(CA) 
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B 



2v 



:ErJ [-< >< > (< U^U" > - <UP ><U^ >) 



(as) 



The partial cancellation leaves: 



A + B 



2v 



-< u'^ Xu'^ > + —< u'' >< u^u'^u'' > 
2o 



1 
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— -< u'' >< u^u^ >< > 



(ae) 



We now evaluate each of the terms in this expression . The first term is 
given by: 



Erj{< U" >< >) = ErJ 



-(x-rV)2/25 



(a?) 

Using the technique of Appendix B, one finds after some work that the only 
non vanishing contribution comes from the a — (5 term in (C.7) and gives: 



1 /Tl 

Er,f{< >< >) = 6 N . 



m 



(as) 



The same technique can be applied to each term in (C.6). One obtains: 

1 — m 



Er,f{< >< uVuP >) = 



S' {N' + 2N) 



iC.9) 



and: 



1 — fji 

Erj{< u" >< u^u" >< u" >) = —3- {N^ + (2 - m)N) {C.IO) 



The sum gives in the end: 



A + B = -- 



1 /2iyV l-m 



5 N 



m 



(au) 



One should notice that the terms of order N"^ vanish automatically due to 
the structure of the velocity field (4.1), and independently from the value of 
m. 

We now proceed to the evaluation of the last term, C, in the correlation 
identity (C.l). This point needs a little care because one must be more 
precise about the correct prescription of the forcing term, whether it is of Ito 
or Stratanovitch type. To settle this issue in a pedestrian but safe way, we 
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have discretized the time in units of tq <^ 1. The Hopf Cole mapping can be 
carried out in this case, and we deduce from this computation: 



C^-f^^f^^-Ne (C.12) 

Using the value (3.12) of 6: 6 = 2mA^, the relation (3.7) between m and 
the Reynolds number: (1 — ^)i-^/2 = ^m^Re^, together with the definition 
(2.4) of Re: Re^ = eA^/i/^ we finally get: 



Re 



\im A+B + C^O (e fixed) (CIS) 



which establishes the correlation identity for the energy balance (C.l). 

Using the same techniques, we have also checked that the correlation 

identities corresponding to ^V.-tT = and ^(('w)^)^ = also hold for any N 
in the limit of large Re. 



30 



Appendix D: Replica computation of the velocity two points 
correlation 



The aim of this Appendix is to compute the two point correlation 



V(^x,t)v{y,t) directly from the replica method of sect. 3. This will provide 
a check that the intermediate physical representation (B.l) is safe. In the 
course of this computation we shall also establish useful replica identities 
which can be of wider interest. To keep the computations simple we consider 
only the one dimensional case, = 1. We start from the result (3.8) for 
the partition function for n replicas of the polymer arriving at time t at the 
points Xa, a = 1, ...,n: 



Z{xi,t)...Z{Xn,t) =C^exp I - ^ Rn{a)n{b)Xa{t)Xb{t) 



{D.l) 



We shall compute Z{x,t)^^'^Z{y,t)'^/'^ which can be written as: 



Z{x,t)^/^Z{y,t)^/^ = cJ2 



a,b 



1- cr„ 



-X + 



-y 



1 - cTfe 1 + o-fe ■ 
— r — X H — y 



[D.2) 



where the Z){o-} is over n Ising spins cr^ = ±1, with the constraint that 
Ea c^a = 0. 

As an intermediate step to this computation, we first compute the action 
of J2'{cr} oiito a polynomial in the spins. Let ai,...,ak be k different replica 
indices. We define: 



{D.3) 



Clearly Ai = 0. As for A2, one can use = J2a,b ^'{a} ^a'^b = n+n(n— 1)^2 to 
deduce A2 = 1/(1 — n). The general result can be deduced from an iteration 
of the above procedure: 



A, 



Ak = 
1 3 



if k is odd 

k-1 



1 — nS ~ n k — 1 — n 
and this can be written in the form: 



if k is even 



in+l 



A, 



dy {cosh{y)Y /2{tanh{y)f 



{DA) 



{D.5) 



B{-n/2, -n/2) J-c 
This expression is well defined for n < 0, and can be continued analytically 



to positive n. B is Euler beta function From this expression of At one 
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can deduce the generating function: 



g[hi...hn\ =X!exp I X! ) = (n.cosh{ha)). 

{a} V a / a 

n 

tanh{hai)...tanh{ha^) 

k=0 ai<...<ak 



•>n+l 



B{-n/2,-n/2) 



dy Y[ cosh{y + ha) 



{D.6) 



a=l 



Formula (D.6) is a replica identity which may turn out useful in other con- 
texts. Here we will use it to compute v{x)v{y). Starting from (D.2), and 
using the fact that J2b ^ab — (derived in Appendix A), we have: 



Z{x, t)-/^Z{y, = c exp ( I 5^ R^baaab 



a.b 



where a is defined as: 



{x - yf 



{D.7) 



(D.S) 



We call as usual r, ri, ro the various elements of the hierarchical Rab matrix. 
Using the fact that ro = 0, we get from (D.6): 



Z(x,i)"/2Z(y,i)"/2 

/oo 
dho 
-oo 



2n+l 



cxp (^a^(r-ri)^ 



- n/m 



B{-n/2,-n/2) 
Dhcosh"'{ho + h^Qxl) 

exp ^ct^(^ — J^i + mri)^ 



5(-n/2,-n/2) 

d/loe"''° J Dh(l + e-2''0-2/»V^)-2mari^ 



n/m 



{D.9) 



where Dh = dh/ y/2n exp {—h^ / 2) . Changing variables to ^ = exp(— 2/?.o) and 
h = m^/fi{z — a;)sign(a; — y), and using the fact that f — ri + mri = 0, we 
derive: 

2n 



Z{x,t)^/'^Z{y,t)^/'^ 



B{-n/2,-n/2) Jo pi 

. — -| n/m 



{D.IO) 



Having computed Z{x,t)'^/'^Z{y,t)'^/'^, the velocity correlation is easily de- 
duced as: 



Z'{x) Z'(y) 



161/2 
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It is straightfoward to check that this rephca computation agrees with the 
direct physical space computation of Appendix B. More precisely we recover 
the expression (B.5),(B.10),(B.ll) for the correlation, using the fact that 
mri6 = 1. 

Figure Captions 

Fig. 1: Typical snapshots of the velocity field as given by Eq. (4.1), with 
m = = 10~^, in one dimension N = 1 (Fig. 1-a), or two dimensions 
N = 2 (Fig 1-b, where we have plotted in grey levels the modulus of the 
velocity field). 

Fig 2: Comparison of the numerically determined velocity correlation 
function (directly from Eq. (4.1)), and our analytical formula, obtained after 
integrating Eq. (4.7), multiplied by m^. Note the linear regime at small x — y. 

Fig 3. Sketch of the 'multifractal' spectrum ({q), giving the r dependence 
of the q th moment of the velocity field, both for the forced Burgers' equation 
(triangles) and hydrodynamical turbulence (circles - see, e.g. |T^). Note that 
C(3) = 1 for both models (Kolmogorov's scaling). Intermittency corrections 
(i.e. the departure of ({q) from |, as given by the dashed line) are much 
stronger in Burgers' turbulence: this is due to the fact that singularities are 
concentrated on A^ — 1 dimensional structures, rather than on vorticity tubes. 
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